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ABSTRACT 


The total-power output from a nuclear reactor can be increased by tailoring the 
axial-power shape if reactor power is limited by fuel- centerline temperature. The po- 
tential power gains and the axial- power shape required to obtain a constant fuel- 
centerline temperature are calculated using steady- state heat-transfer calculations on a 
long fuel rod cooled externally by forced convection. Steady-state two-group neutron dif- 
fusion theory is used to predict the axial-fuel loading required to obtain the desired axial- 
power shape. Fuel loadings and neutron-flux distributions are presented for various 
axial-power shapes required in a reflected, water moderated, and cooled reactor. 
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AXIAL POWER TAILORING TO OBTAIN CONSTANT FUEL-CENTERLINE 
TEMPERATURE IN A NUCLEAR REACTOR 
by Harry W. Davison 
Lewis Research Center 

SUMMARY 

When the power output from a nuclear reactor is limited by the temperature of the 
fuel, the power can be increased by tailoring the axial-power distribution to produce a 
constant fuel- centerline temperature. Heat-transfer and neutronic calculations have 
been made to determine (1) the potential power gain which can be obtained by power tail- 
oring, (2) the axial-power shape required, and (3) the axial-fuel distribution required to 
obtain the required axial-power shape. 

The steady-state heat-transfer calculations were made for a long fuel rod cooled by 
forced convection to an incompressible coolant. The neutronics calculations are based 
on the one-dimensional diffusion theory using two energy groups. The results of this 
study indicate that the power of an unreflected, axially unzoned reactor, limited by fuel- 
centerline temperature, could be increased by at least 50 percent by tailoring the axial- 
power shape. When the convective film coefficient is constant, the required axial-power 
shape is exponential. An exponential-power shape can be obtained by appropriate design 
of the end reflectors and by adjusting the axial-fuel distributions. When the neutron re- 
flectors at the coolant inlet and outlet ends of the reactor are identical, the maximum- 
fuel loading is required at the coolant inlet end of the reactor. The maximum and aver- 
age fuel concentrations required to maintain reactor criticality are larger for small 
reactors than for large reactors . 

Although the required axial-power shape and fuel distribution were calculated for a 
system having a constant film coefficient, the change in the required axial-power shape 
caused by small axial variations in the film coefficient was estimated. 


INTRODUCTION 

In power reactors it is desirable to extract the maximum power per unit volume of the 
system for most efficient operation. The thermal-power capability of a reactor depends 



on the magnitude of the operating limits. These limits are imposed on parameters such 
as fuel or component temperatures, coolant temperature, coolant flow rate, or heat flux 
to permit safe operation of the reactor. If reactor power is limited by fuel temperature 
or heat flux, which depend on the power shape in the reactor, reactor power can some- 
times be increased by tailoring or shaping the power distribution. 

One method of power tailoring relies on distributing or zoning the nuclear fuel within 
the reactor to obtain the desired power shape. Fuel distributions (fuel zoning) required 
to produce a flat radial power distribution have been investigated by Goertzel and Loeb 
(ref. 1) and by Bussard and Delauer (ref. 2). Barth and Haling (ref. 3) investigated 
axial-fuel distributions required to obtain an optimized-power distribution based on crit- 
ical heat-flux requirements in a boiling water reactor. 

In this report, the focus of attention is on axial-power distributions required when 
reactor power is limited by fuel- centerline temperature. If reactor power in an unzoned 
reactor is limited by fuel- centerline temperature, there is usually only one spot in a fuel 
rod where the centerline temperature reaches the temperature limit. At all other axial 
positions along the fuel centerline, the temperature is below the temperature limit. The 
total fuel-rod power and reactor power can be increased by adjusting the axial-power 
shape so that the fuel- centerline temperature reaches the temperature limit and is con- 
stant along the entire length of the fuel rod. Axial-power tailoring to obtain a constant 
fuel-centerline temperature is desirable to obtain higher reactor powers. Axial-power 
tailoring would also allow a reduction in reactor size while maintaining the same reactor 
power, providing the reactor is not criticality limited. 

The purpose of this report is 

(1) to illustrate the potential gains in reactor power that might be obtained by tail- 

oring the axial-power distribution when reactor power is limited by fuel- 

centerline temperature 

(2) to determine the axial-power distribution (hereinafter called the optimum axial- 

power shape) required to produce a constant fuel-centerline temperature 

(3) to determine the fuel distribution required to produce the optimum axial-power 

shape. 

The potential power gains, as a result of axial-power shaping, and the optimum 
axial-power shape are determined as a result of an analytical, steady-state, heat- 
transfer analysis on a fuel rod which is cooled externally by forced convection to an in- 
compressible coolant. The fuel-material distribution required to produce a constant 
fuel- centerline temperature is calculated using the one -dimensional, two-group, steady- 
state, neutron diffusion theory. 


2 



HEAT-TRANSFER ANALYSIS 


K 

/ 


The total reactor power that can be obtained with and without a tailored axial-power 
distribution and the axial-power distribution required to produce a constant fuel- 
centerline temperature are derived in appendix A. The potential gain in reactor power is 
based on the following reactor model: 

(1) The fuel element is a long cylindrical rod which is cooled externally by incom- 
pressible, coaxial, forced- convection flow. The fuel element operates at steady state 
and is heated internally by fission of the fuel material. The convective film coefficient 
and thermal conductivity of the fuel are independent of axial position. 

(2) The radial power distribution in the reactor is constant and is unaffected by axial- 
power tailoring. 

(3) Reactor power is limited by fuel- centerline temperature, that is, the maximum 
reactor power is obtained when the temperature of the hottest spot in the hottest fuel 
element in the reactor reaches a prescribed limit. 

The total reactor power that can be obtained with a tailored axial-power shape is 
given in equation (A31) of appendix A. 
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(A 31) 


The dimensionless parameter 0 is defined in equation (A13). It is a measure of the 
sum of the temperature rise across the convective coolant film and across the fuel rod 
relative to the total coolant-temperature rise across the reactor. Beta depends on the 
coolant flow rate and heat capacity, the convective film coefficient, the radial power dis- 
tribution in the fuel rod, and the dimensions and thermal conductivity of the fuel. In this 
study, 0 will be treated as a thermal-hydraulic-design parameter which influences the 
axial-power distribution, and, consequently, affects total reactor power. As indicated 
in equation (A13), 0 could vary from zero to infinity. For most pressurized water 
reactors, for example, 0 is about 10. 

The reactor power obtained with a constant fuel distribution (untailored axial-power 
shape) is given by equation (A39). 
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where is a function of the parameters 0 and the axial-power shape; that is, the 
values of Xj, Xg, and Xg depend on the axial-power shape. The axial-power shape in 
an unreflected reactor having a uniform fuel distribution is described by a cosine function 
(xj = 0 and Xg = x^ = 1). When neutron reflectors are added to both the coolant inlet and 
outlet ends of the reactor, however, the axial-power shape is flatter and can be approxi- 
mated by a "chopped cosine" (x^ > 0 and Xg < Xg). The chopped cosine axial-power 
shape is normally described by the ratio of the maximum to average power generation 
rates n^ rather than in terms of x^, Xg, and Xg. 

With the same fuel-temperature limit, coolant flow, and coolant outlet temperature, 
the power that can be obtained by power tailoring relative to that which would be obtained 
without power tailoring is 



Because the coolant flow and outlet temperature are the same, the power gain is achieved 
by reducing the coolant-inlet temperature. 

The ratio P e /P^ is shown in figure 1 as a function of 0 and n . The potential 
power gain that can be achieved by power tailoring an unreflected reactor having a con- 
stant fuel distribution is described by the curve labeled il = 1. 57. The other curves 
(np < 1. 57) represent the potential power gain that can be achieved by power tailoring 
reflected reactors initially having constant fuel distributions and a ratio of maximum-to- 
average power generation rates of np. The potential power gains are less for large 
values of 0 than for small values of 0. However, even for large values of 0 (0 > 10) 
tailoring the power shape in an unreflected reactor would allow about a 50-percent in- 
crease in reactor power. When 0 = 10 (as it is in many pressurized water reactors) or 
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Figure 1. - Potential gain in reactor power which might be achieved by axial power tailoring when reactor 
power is limited by fuel centerline temperature. Reactor powers are calculated for same coolant flow, 
coolant effluent temperature, and maximum fuel centerline temperature. 


greater, the potential power gain obtained by axial-power tailoring is less for reflected 
reactors than for unreflected reactors. For example, the power of an unreflected 
(np = 1.57), pressurized water-moderated, reactor, which is limited by fuel- centerline 
temperature, could be increased by at least 50 percent by tailoring the axial-power shape 
to produce a constant fuel- centerline temperature. If this reactor had neutron reflectors 
at the coolant inlet and outlet ends such that n p = 1. 10, the power could be increased by 
only about 10 percent by axial-power tailoring. The reactor power gain is smaller in the 
second example because the initial (untailored) reactor power of a reflected reactor is 
higher than the initial reactor power of an unreflected reactor. The final (tailored) 
reactor power would be the same in both cases . 

These indicated large power gains may not actually be achievable because reactor 
power may be limited by some parameter other than fuel temperature before the calcu- 
lated power gain can be achieved. Although the power gained by axial-power tailoring 
may not be as large as that indicated in figure 1, a lesser increase in power would still 
be an advantage. 

Having determined the potential reactor-power gains that might be achieved by axial- 
power tailoring, what is the axial-power distribution required to produce a constant fuel- 
centerline temperature? The required axial-power shape is exponential, and is defined 
in equation (A29). 
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This power shape is shown in figure 2 for various values of 0. The maximum power is 
required at the coolant- inlet end of the reactor where the coolant temperature is lowest. 
The minimum power is required at the outlet end where the coolant temperature is 



Figure 2. - Axial -power distribution required for constant 
fuel-temperature distribution. 


highest. The power shape required for larger values of £ is flatter than that required 
for smaller 0 values. This is because larger values of 0 indicate smaller coolant- 
temperature rises relative to the temperature rise across the fuel and convective coolant 
film. Therefore, only a slight reduction in the axial power from inlet to outlet is re- 
quired to compensate for the small coolant-temperature rise. 

The previous discussion is based on the assumption that the parameter 0 is inde- 
pendent of axial position. Actually, the film coefficient may vary along the length of the 
fuel element as the coolant temperature increases. The axial-power distribution re- 
quired for small axial variations in film coefficient is estimated with equation (A24) of 
appendix A. 
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The factor a ^ is a measure of the total change in the film coefficient along the length of 
the fuel element, and it can be calculated using equation (A22). The change in the re- 
quired axial- power distribution caused by small variations (cuj « 0) in the film coeffi- 
cient is illustrated in figure 3 where the ratio Y(x,/3,aj)/Y(x,/3, 0) is presented as a 



Figure 3. - Effect of total change in film coefficient along 
length of fuel rod on required axial-power distribution. 


function of axial position and For many reactors a j is small (-0. 1 < a j < 0. 1). 

Positive values of a ^ result when the convective film coefficient decreases with tem- 
perature. This means that the heat-transfer resistance is lower at the coolant- inlet end 
than at the outlet. Therefore, the power at the inlet must be higher and that at the outlet 
lower than would be required with a constant film coefficient (a j = 0) . The result is a 
steeper power shape when > 0. Conversely, negative values of a j are caused when 
the film coefficient increases with temperature. The required axial-power shape is 
flatter than that required with a constant film coefficient because a smaller reduction in 







power is required to compensate for the coolant-temperature rise. 

The previous discussion has illustrated the potential reactor power gains that might 
be achieved by tailoring the axial-power shape and has established the axial-power shape 
needed to obtain a constant fuel- centerline temperature. The next problem is to deter- 
mine whether the axial-power shape can be tailored to produce a constant fuel- centerline 
temperature and how this tailoring can be done. In this analysis it will be assumed that 
the film coefficient is constant. When the film coefficient is constant, the required 
axial-power shape is an exponential as indicated in equation (A28). 


NEUTRONICS ANALYSIS 


There are several methods of power tailoring available to the reactor designer. 

The most common methods are achieved by design of a neutron reflector or by appro- 
priate distribution of fuel or absorber materials. 

In this study, power tailoring will be attempted by adjusting the fuel distribution 
(fuel zoning) and by reflector design. 

Physically, the fuel distribution could be adjusted by varying the isotopic proportion 
of fissionable material in the fuel rod. For example, if uranium is used as fuel, the 
proportion of uranium 235 to uranium 238 could be adjusted to obtain the required dis- 
tribution of fissionable material in a constant-diameter fuel rod. hi the following dis- 
cussion the required fuel distribution is expressed in terms of the macroscopic fission 
cross section S^. The fuel distribution required to obtain an exponential axial-power 
shape in a reflected reactor is derived in appendix B. 

The neutronics calculations are based on a one-dimensional steady-state diffusion 
calculation. Two energy groups of neutrons are assumed, such that all fissions result 
from the absorption of thermal neutrons. The fission neutrons all appear with the same 
energy, which is higher than the thermal energy. These fast neutrons are lost from a 
unit volume of reactor by diffusion and by slowing down to thermal energies. The 
thermal neutrons are formed by the slowing down of the fast neutrons. They are lost 
from a unit volume of reactor by diffusion and absorption in both fuel and nonfuel mate- 
rials. The effect of neutron reflectors located at both ends of the reactor is simulated 
by specifying neutron albedos for both the fast and thermal fluxes. The fast and thermal 
neutron fluxes are represented by equations (B14) and (B22) in appendix B. 
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The parameters y^ and y can be related to the heat-transfer parameters through 
equations (B37) and (B35). 
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0(X 2 - Xj) 

Jo rH(r)dr 
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V = [iS(X 2 - Xj)] 


(B35) 


The coefficients Aj, A 2 , Ag, and A^ are given in equations (B17), (B18), (B27), 
and (B28). The coefficients Aj and A 2 depend on the core length, the fast-neutron 
properties of the core, and the fast-neutron albedo at both ends of the reactor. The 
coefficients Ag and A^ depend on Aj and A 2 , the thermal properties of the core, 
and the thermal-neutron albedo at both ends of the core. The fuel distribution required 
to produce an exponential axial-power distribution is 


S f (X) = 


D 1 y 1 exp(-yX) 

p<P 2 (*) 


The required fuel distribution depends on 

(1) the fast and thermal- neutron diffusion coefficients Dj and D 2 

(2) the thermal-neutron absorption cross sections of nonfuel materials 

(3) the removal cross section for fast neutrons S R 
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(4) the number of neutrons produced as a result of thermal-neutron absorptions in 

the fuel r] 

(5) the length of the reactor Xg - Xj 

(6) the fast and thermal neutron albedo at both ends of the reactor Kj j, Kg j, 

K l,2> K 2, 2 

(7) the exponential coefficient used to describe the axial-power shape y. 

Because of the large number of independent variables involved (eleven) in the deter- 
mination of fuel tailoring, it is desirable to reduce the number of variables by selecting 
a specific type of reactor for analysis. A water moderated and cooled reactor was 
selected for this analysis. The reactor is fueled with uranium 235, and contains alumi- 
num structural components. The following reactor parameters, obtained from refer- 
ences 4 and 5, were used in this study: 


Fast neutron diffusion coefficient, Dj, cm 1.42 

Thermal neutron diffusion coefficient, Dg, cm 0.262 

Fast neutron diffusion length, L^, cm 5.95 

Thermal neutron diffusion length, Lg, cm 4.10 

Neutron removal cross section, cm - * 0.0398 

Number of neutrons produced per thermal neutron absorbed in fuel, t] 2.08 

Number of neutrons produced per fission, v 2.44 



Figured - Neutron flux distribution in reac- 
tor with exponential axial power shape 
Y(x) = Yjexp -(x/p). Reactor length, 

120 centimeters; no thermal leakage from 
reactor {Kj 2 = K 2 2 = 



Figure 5. - Fuel distribution required to 
produce exponential axial -power distribu- 
tion Y<x) = Yjexp-fx/p). Reactor length, 
120 centimeters; no thermal leakage from 
reactor 2 = K 2 2 = 
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It will be assumed that the reactor length is 120 centimeters and that it has a moder- 
ating reflector at both ends so that the fast-neutron albedos are zero (K^ j = Kg j = 0). 

The axial neutron-flux distribution and the axial-fuel distribution required to main- 
tain a constant fuel- centerline temperature are illustrated in figures 4 and 5 for three 
different axial-power shapes (3 values of /3). The fast and thermal neutron-flux distri- 
butions for three different axial-power distributions are shown in figure 4. The re- 
flector at the coolant- inlet end of the reactor and the reflector at the outlet are identical 
and are designed such that there is no thermal-neutron leakage from the reactor 
(K^ g = K 2 2 = *)• -Although the two reflectors are identical, the neutron-flux distribu- 
tions are asymmetric. The asymmetry in the flux distribution is caused by the asym- 
metry in the fuel distribution (see fig. 5) required to produce an exponential- power shape. 
More fuel must be added to the inlet end than to the outlet end because a higher power is 
required at the inlet end. The fuel distribution required to produce a flat power shape is 
symmetric and is lower in the center than at either end. Although in this example the 
fuel content remains constant over the central 60 percent of the length of the fuel rod, 
this region of constant-fuel content depends on the length of the fuel, the neutron removal 
cross section, the diffusion coefficient, and the diffusion length. 



Figure 6. - Fuel content required to produce expo- 
nential axial -power distribution Y(x) - Yj 
exp-(x/j}). No thermal leakage from reactor 
(K| 2 ” ^2 2 ~ 
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It is desirable to learn what fuel distributions are required for different reactor 
lengths. The maximum and axially averaged values of fuel content for various reactor 
lengths are illustrated in figure 6. Again the thermal albedos at the inlet and outlet ends 
of the core are one. The smallest reactors that can be made critical while maintaining 
an exponential axial -power shape (Y(x) = Yj exp(-x/jS) are represented by the vertical 
lines. For example, the smallest reactor that could be made critical while maintaining 
an axial-power shape of the form e“ x (/3 = 1) is about 27 centimeters long. 

Typically, smaller reactors require larger fuel loadings than larger reactors. For 
example, the maximum-fuel content and the axially averaged fuel content required to 
obtain an axial-power shape of the form e” x (/3 = 1) in a 50- centimeter- long reactor are 
0. 11 and 0. 023 reciprocal centimeter, respectively. The maximum and average fuel con- 
tent required to maintain the same axial-power shape in a 100-centimeter-long reactor 
are 0. 078. and 0. 017 reciprocal centimeter, respectively. 

If a reactor is to be designed with an exponential- power shape and if the maximum 
fuel content is limited, the reactor length can be increased to obtain the desired 

power shape. If, however, an exponential-power distribution is desired in an existing 
(fixed length) reactor, it may be possible to increase the fuel content to obtain this 

power shape. 

The previous discussions assume that there is no thermal-neutron leakage from the 
reactor (K^ 2 = K 2 2 = T ^ e e ^ ec ^ °* the reflectors on the maximum and average fuel 
content required are illustrated in figures 7 and 8, respectively. The required fuel con- 
tent is given as a function of the thermal-neutron albedo at the inlet Kj g an< * outlet 
K 2 2 enc * s °f the reactor. When the reactor is unreflected or when the two reflectors 
are identical, the maximum-fuel content is required at the inlet ends of the reactor. As 
more reflector is added to the inlet end of the reactor (Kj 2 increases), a point is 
reached where the maximum-fuel content will be required at the outlet end. For example, 
in figure 7 for no reflectors (Kj g = Kg 2 = ^)> max i mum -fuel content (about 0. 53 
cm"*) is required at the inlet end of the reactor. As more reflector is added to the inlet 
end, Kj g can increased to 1.0, and there is no thermal-neutron leakage from the 
inlet end. The fuel content required at the inlet is about 0. 07 reciprocal centimeter. 
However, the fuel content required at the outlet end is about 0. 33 reciprocal centimeter 
because no reflector has been added there. Further improvements in the inlet- end re- 
flectors have little effect on the fuel content required at the outlet end. The outlet-end 
reflector must be modified to increase Kg g before the maximum-fuel loading can be 
reduced further. It should be noted that the thermal- neutron albedo is limited: it de- 
pends on the reflector size and material, but can generally have values between zero and 
slightly greater than one. 


12 




Figure 7. - Effect of reflectors on maximum fuel con- 
tent required to produce exponential axial power 
distribution Ylx) = Yj exp - x. Reactor length, 

120 centimeters; thermal -hydraulic factor, (3, is 1. 



Figure 8. - Effect of reflectors on average fuel content 
required to produce an exponential axial power dis- 
tribution Ylx) = Yjexp - x. Reactor length, 

120 centimeters; thermal -hydraulic factor, 1. 
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CONCLUSIONS 


When reactor power is limited by fuel- centerline temperature, the power can be 
increased by tailoring the power distribution. The power gain depends on the operating 
conditions. However, the power in an unreflected reactor could be increased by at least 
50 percent if the optimum axial-power shape can be obtained. The optimum axial-power 
shape yields a constant fuel-centerline temperature. When the convective -film coeffi- 
cient in the reactor is constant, the optimum axial-power shape is of exponential form. 

Generally, the axial-power shape required to maintain a constant fuel- centerline 
temperature is flatter when the coolant-temperature rise is small. 

The axial-power distribution can be tailored by proper distribution of the fuel and by 
proper design of reflectors at the coolant inlet and outlet ends of the reactor. When the 
two reflectors are identical, the maximum-fuel loading is required at the coolant inlet 
end of the reactor. As with unzoned reactors, larger fuel concentrations are required 
for small reactors than for large reactors . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 10, 1968, 

120-27-06-17-22. 
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APPENDIX A 


■ 


AXIAL POWER DISTRIBUTION REQUIRED TO OBTAIN 
CONSTANT-FUEL-TEMPERATURE DISTRIBUTION 

The axial-power distribution required to produce a constant fuel-centerline tempera- 
ture is determined by calculating the axial-temperature distribution as a function of the 
unknown axial-power distribution. The functional form of the temperature distribution is 
differentiated with respect to the axial distance; this derivative is allowed to vanish (con- 
stant temperature), and the resulting differential equation is solved to determine the 
optimum axial-power distribution. 

The axial-temperature distribution is calculated for a reactor fuel rod (see fig. 9) 
operating at steady state and cooled by an incompressible coolant. The fuel assembly is 
cylindrical, the thermal conductivity is constant, and the length is much greater than the 
diameter. The temperature distribution in the fuel element is described by the Poisson 
equation: 



Figure 9. - Reactor fuel rod. 
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d 2 e + 1 de + _ q(r,x) 

SR 2 R 0R 3X 2 k M 

This equation can be nondimens ionalized to unit order by the following substitutions: 


a 

max 



AX 


The governing equation in terms of the dimensionless parameters is 

®max 0 ^ + ^max ^ + 0 max 0 ^ _ Q(R,X) 

R 2 01,2 R o r o 0r ( AX ) 2 0x2 k]VI 

S 2 T + 1 ST + / R o\ 2 8 2 T = R o q < r > x > 

9r 2 r 0r \ AX / 3x 2 0 max k M 

Where q(r,x) is Q(R,X) expressed in terms of the nondimensional variables r and x. 
All variables are of unit order, and R q /AX « 1. When the axial-heat loss is small 
relative to the radial-heat loss, equation (A2) can be approximated by 

8 2 T + 1 ST = R p q ( r > x ) 

Sr 2 r 0r 0 max^M 

For a solid rod, heated internally and cooled externally, the following boundary condi- 
tions are used: 


^1=0 at r = 0 
dr 
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T = T(l,x) 


at r = 1 


The temperature anywhere within the rod is 


T(r,x) = T(l,x) + 


— f/'fx 

k M 0 max (-4 


r rq(r,x)drj^ r ^ 


The temperature along the axis of the fuel rod is 


T(0,x) = T(l,x) +■ 


-S-f r\r 

k M 0 max L-'O 


rq(r,x)dr — 

J r 


(A4) 


(A 5) 


Let the volumetric heating rate be expressed as the product of a function of r,H(r) 
and a function of x, Y(x). This relation is substituted into equation (A 5). 

RqY(x) 

T(0,x) = T(l,x) + — 

k M^max 



The temperature at the surface of the fuel rod T(l,x) is equal to the local bulk- 
coolant temperature T £ (x) plus the temperature rise across the convective coolant film. 
The local bulk-coolant temperature can be expressed in terms of the coolant-temperature 
rise. 


T c (x) =T 1 + m(x)(T 2 - Tj) 

where m(x) is the fraction of the total fuel-rod power generated up to axial position x. 


m(x) = 


f X Y(x)dx f l rH(r)dr f* Y(x)dx 

yxj Jo *'x 1 

f* 2 Y(x)dx f 1 rH(r)dr f* 2 Y(x)dx 

yXj •'0 ✓Xj 
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The temperature rise across the convective coolant film can be expressed by 


T(l,x) - T (x) =(Q \ 

\A/h0 max 


where n(x) is the ratio of the local to axially averaged power generation in the fuel rod. 


n(x) = 

H(r) f 2 Y(x)dx 


Y(x) 


\ f X2 Y(x)dx f X 

yxj Jx j 


Y(x)dx 


Therefore, the temperature at the surface of the fuel rod T(l,x) can be expressed by 


T(l,x) =Tj + 


f 

yxj 


Y(x)dx 


(T 2 - Tj) + 


FC p (T 2 - T 1 )Y(x) 


(A7) 


f * 2 Y(x)dx 2irR hAX f 2 Y(x)dx 

Jx 1 Jx 1 


When equation (A7) is substituted into equation (A6) 


T(0,x) - T 


f 

1 Jx l 


Y(x)dx 


FC p Y(x) 




r rH(r)drl— 


t 2 - t 


1 f 
JX 1 


2 Y(x)dx 2 ttR 


o-x r 

Jx i 


Y(x)dx 


k M^ T 2 " T l^ 0 


max 


(A8) 


The total coolant-temperature rise can be expressed by 


T 2 " Ti = 


27tR 2 AX /* Xq /*1 

, = 2 / 2 Y(x)dx / rH(r)dr 

1 FC_0 „ Jx * JO 

p max 1 


(A9) 


When equation (A9) is substituted into equation (A8), the result is 
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r 

T(0,x) - Tj yxj 


r 


Y(x)dx 


FC p Y(x) 


t 2 -t 


1 f* 2 Y(x)dx AX f* 2 Y(x)dx 

Jx J */Xj 


, j\fo rH(r)dr ]~ 

2lrk M f 1 rH ( r ) dr 
J 0 


27rR Q h 


and 


r 

T(0,x) - Tj Jx 1 


Y(x)dx 


FC p Y(x) 




t 2 -t 


<- + 
h 


1 f* 2 Y(x)dx 2irR AX f * 2 Y(x)dx h k Mjf 1 rH ( r ) dr 

yxj y Xj 


(A 10) 


This equation can also be expressed in terms of dimensionless Reynolds, Prandtl 
and Nusselt numbers. 


x) T r Y<x>dx f 

> x ' ~ T 1 _ yx l + RePrSDeY(x) J _1_ + JQ 

1 f * 2 Y(x)dx 8irR AX f * 2 Y(x)dx ^ k M De rH ( r ) dr 

Jx x Jx 1 l 


t 2 -t 


(All) 


Equation (All) can be rewritten in the following form: 


/ X Y(x)dx 

. A ^ 


t 2 -t 


+ /3 


Y(x) 


(A 12) 


1 / 2 Y(x)dx f 2 Y(x)dx 

•OC-. ./x. 
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where 


FC 


/3 = 


'•/‘U 


2ttR 0 AX 


r rH(r)drJ— 'j 


— + 
h 


“m / 1 rH < r > dr 


(A 13) 


/3 


RePrSDe 

87tR q AX 




r rH(r)drj— 


Nu 


k M De^T 1 rH(r)dr 


(A 14) 


The axial-power distribution Y(x) required for a constant fuel- centerline tempera- 
ture is determined by allowing the derivative of the temperature to vanish. 


dT (°’ x > = 0 = Y(x) (l + M +13^) 
dx \ dx/ dx 


(A 15) 


The coolant properties depend on the coolant temperature. Therefore, the film coeffi- 
cient and 0 depend on the coolant temperature. 


dg = M (de\ 

dx \d0 / \dx / 


(A16) 


The bulk coolant temperature at any axial position can be expressed by 


2ffR?AX 


e = e 1 + 


— f 1 rH(r)dr f x Y(x)dx 

FC d *^ X 1 


(A17) 


d 9 
dx 


2nR^AX 


f rH(r)dr 

FCL JO 


Y(x) 


(A 18) 


The functional relation between 0 and temperature can be approximated by 
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(1) calculating values of /3 for various values of coolant temperature 

(2) fitting the values of /3 to a polynominal in coolant temperature. 

/3 = a Q + aj0 + a 2 0 2 + . . . + a n 0 n (A19) 

= a.-, + 2a 0 9 + 3a« 0 2 + .. . + naj"' 1 (A20) 

d0 1 2 6 n 

When equations (A 18) to (A20) are substituted into equation (A15), the result is an 
equation describing the axial-power distribution Y(x) as a function of distance along the 
coolant channel x in integral-differential form. This equation can be solved in closed 
form when j3 is constant. For most incompressible coolants such as water and liquid 
metals, (3 can be approximated, over the temperature range of interest, by either a 
constant or a linear function of coolant temperature. When j3 is a linear function of 
coolant temperature, 


M 

d<9 


= a l 


When equation (A18) to (A20) are substituted into equation (A15), the differential equation 
describing the axial-power shape is 

/3 ^ + [l + 0! jyJ Y - 0 (A21) 


where 


“1 = 


2irajR^AX 


FC 


r 1 

Jo 


rH(r)dr 


(A22) 


The meaning of aj can be more easily recognized by substituting the temperature de- 
rivative of equation (A13) into equation (A22). The temperature derivative of equa- 
tion (A 13) is 

«£- ai - - FC p 

de 1 27rR 0 AX dT \h/ 
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When this is substituted into equation (A22) , the result is 


°* - M 1 rH H« © 

Therefore a j is a measure of the change in the film coefficient along the length of the 
fuel rod. Comparing equations (A15) and (A21) results in 

^ = a'-.Yfx) 
dx 1 

dp = a-. f x Y(x)dx 

\ Vx x 

The function /3 can then be expressed as 


0 = 0i 


+ a -i f x Y(x)dx 


where /3j is evaluated from equations (A 13) or (A 19) using physical property data eval- 
uated at the coolant- inlet temperature. When the axially averaged power distribution is 
normalized to one 


P Y(x)dx 

— •'Xi 

Y = t = 1 

x 2' x l 


By definition, 

x 2 - x i = 1 

Therefore, 


f* 2 Y(x)dx = 1 
•/x.. 
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0 < f x Y(x)dx < 1 
Jx l 


When o!^ « /3-p the solution of equation (A21) can be approximated by 

>Y 


f dY = - f x 

J Y(1 + ajY) fiA, 


dx 


the solution is 


Y(x) = Yj_ 


(1 + o; 1 Y 1 )exp ^ t ^ - o-jY 


r 


where Yj is the relative power at the coolant- inlet end of the reactor. Since 

f * 2 Y(x)dx = 1 

Jx l 

Yj can be expressed in terms of and /3 as follows: 

1 - expCi) 

Y. = W. 


Y(x) _ 


“{“’(■ «j’ 1 ] 

expQ - 1 


Y 

(a 

A A 

"1 

1 - expl — 

\/3 

r) - exp (- 


(A23) 


(A24) 


The power distribution in the fuel rod that yields a constant fuel- centerline temperature 



q(x,r)=Y 1 H(r 


{ 


(1 + a-j^Y^exp 




1-1 


“ 1 Y 1 


With this power distribution, the coolant-temperature rise (0 2 " 0j) and reactor power 

P may be expressed as functions of the fuel temperature limit 0 L using equation (A 12). 
© 



( 

1 \ 


r 

/a A 


( All 

d 2~ d l~( d L~ 9 2* a l 

exp (- - 

' i 

- - 1 

(3/ J 

< 

B 

L 

L 1 - -{ t )\ 

+ a i 

[_ H J 



(A25) 


(A26) 


When the film coefficient does not vary along the length of the fuel 

= 0 
dX 


and 


Y 


1 




(A27) 


The axial power shape required to produce a constant fuel-centerline temperature is of 
the following exponential form. 


(x - x,) 

Y(x) = Y, exp 

1 /3 



(A28) 


(A29) 
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The coolant temperature rise and reactor power are 


9 


2 


P 


e 


" d l = ( 0 L " d 2 > 


— (9 T - 9 n) 


exp 



(A30) 


(A31) 


The reactor power that can be achieved by power tailoring can be compared with the 
power from a reactor having a constant fuel and neutron absorber distribution. A con- 
stant fuel and absorber distribution yields an axial-power shape which is symmetric 
about the reactor midplane and can be described by a cosine function when the origin of 
the coordinate system is located at the reactor midplane. In this study, however, the 
origin of the coordinate system is assumed to be above the coolant- inlet end of the re- 
actor as shown in figure 9. Although the axial-power shape is called a chopped cosine to 
conform with convention, it is described in equation (A32) as a sine function. 


Y(x) = n p 



n 


P 




dY(x) _ 
dx 



(A 32) 


(A33) 


The maximum-fuel temperature and reactor power in a reactor having a chopped 
cosine power shape is calculated using the same methods used to derive equation (A31). 
hi this case, however, the axial-power shape is known and the position of maximum 
temperature x must be determined in order to evaluate the maximum-fuel temperature 
from equation (A 12). When equation (A32) is substituted into equation (A 12) and evaluated 
at the position x of maximum temperature the result is 

hr 
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(A 34) 



The position of maximum temperature is determined by substituting equations (A32) and 
(A33) into equation (A 15). It is assumed that the film coefficient is constant, that is, 
d/3/dx = 0. The result is 



Using this equation and the trignometric identities yield 



(A 3 5) 


(A36) 


The maximum temperature is determined by substituting the values from equa- 
tions (A3 5) and (A36) into equation (A34). 



(A 37) 
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With a cosine power distribution, the coolant-temperature rise (0 g - 0j) and reactor 
power can be expressed as a function of the fuel-temperature limit djj 


I 

i 


e 


2 




P 


d _ 



(A38) 


(A39) 
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APPENDIX B 


AXIAL-POWER TAILORING BY VARYING AXIAL-FUEL DISTRIBUTION 

The method used to determine the fuel distribution required to obtain an exponential 
power shape is similar to the one discussed by Goertzel and Loeb (ref. 1). The axial 
neutron-flux distribution is calculated as a function of the known exponential source 
(power) distribution using a one-dimensional, twp-energy group neutron-diffusion model. 
The required fuel distribution is proportional to the power distribution divided by the flux 
distribution. The reflector is simulated by specifying the neutron albedo at both ends of 
the reactor. The neutron population in a unit volume of the reactor can be approximated 
for each energy group by 


d 2 cpAX) 

Dj 1 Sr^X) + i/S f (X)? 2 (X) = 0 (Bl) 

dX 2 


Dr 


d V 2 (X) 


dX 


- S >2< X > 


s£(X)<P 2 (X) + Sr^X) = 0 


(B2) 


where the subscripts 1 and 2 refer to fast and thermal energy groups, respectively. In 
the fast group, neutrons are produced in a unit volume by fissions and are depleted from 
a unit volume by leakage and by slowing down in the moderator. In the thermal group, 
neutrons are produced in a unit volume by slowing down the fast neutrons. They are de- 
pleted by leakage, absorption in the fuel material and absorptions in the nonfuel mate- 
rials. The neutron diffusion equations can be rewritten as a function of the desired 
source (power) distribution 


d 2< ? j(X) (p^X) 

dX 2 L 2 


-G(X) 


d * 2 « X) ^ fS, l( x)=G l( X) 


dX 


D 9 ?? 


(B3) 


(B4) 


where 
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2 D 1 


G(x) = 


i/S f (x)<P 2 (x) 


r 2 _ D 2 

J-io 


The solution to the first equation is 


(p^X) = <pJ(X) + <pf(X) 


The first term on the right is the solution of the homogeneous equation. The second term 
is a particular solution. The solution of the homogeneous equation is 


*?<X) = A - 1 exp + A 2 


(BIO) 


The particular solution can be found using operator notation. 


° 2 " — K (X) = - G( X ) 
, l 2 i/ 


where 


^2 d 
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Let 




— W?(X) = -G(X) 


u = 







(Bll) 



The lower limits, X Q and X Q0 , may be arbitrarily selected to yield the most convenient 
solution. Any particular solution is acceptable provided it satisfies equation (B3). 
Therefore, 


■ *■ • Ai ~t) * "'(f)/ 

A oo A o 


X 


G(X)exp/— \dX dX 


V 1 


(B12) 


When G(X) is an exponential of the form 
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G(X) = r 1 exp(-rX) 


(B13) 


The following limits are selected: 


Xq = 0 and Xqq = °° 


The solution of equation (B12) is 


cp-,(X) = A-, expf— | + A 9 exp (- — | - 

VV \ L i/ 


Lfvi exp(-vX) 
Ljx* - 1 


where Ag is equal to A^ plus terms derived from the particular solution (X) . 

The boundary conditions are established by specifying the albedo at the coolant- inlet 
end of the reactor j and at the coolant-outlet end Kg j. The neutron albedo for any 
energy group is defined as the ratio of the neutron partial current entering the reactor to 
the neturon partial current leaving the reactor. 

Kj j = (B15) 

1.1 J_(Xj) 

J < X 2> , , 

K 9 , = — (B16) 

2.1 W 

J + (Xj) = Kj jJJXj) 

J -( X 2 )= K 2,1 J + ( X 2) 


The neutron current can be expressed by 


j (x) = - Dl <MX) 

+ 4 2 dX 


CD (Xl ^ 1 ^ ^ X ^ 

J (X) = + 


l 



When these equations are substituted into equations (B15) and (B16), 

,, „ , D i 




*P l(Xo) 

(i-k 2 j) — L_f_ = -(i + k 2 j) 

4 *» 1 2 dX 


2 dX 

Dj d<P x (X 2 ) 


These equations are solved to determine Aj and Ag. 


A 2 " 


^1 L 1 
r 2 L? - 1 


exp(-yXj) -^ 3^4 ex p 

.^2 - X 1 

^2^5 " ^4^6 exp 

■- 2 (X 2 - Xj)' 

L 1 


+ J-X, 




(B17) 


A 1 " 


T 2 

^1 L 1 J 

i^ 3 ^ 5 exp(-yX 2 ) -i/'jtf/g exp 

’ / X 2 - X 1 
_\ L 1 

v 2 L 2 - 1 

^2^5 " ^4^6 exp 

L 

" 2(X 2 - Xjf 

L 1 . 


+ yX, 


exp^ — — ^ (B18) 
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+1 = 

^3 = 
^4 = 

*5 = 

^6 = 


^Kl.l 

, d+K 1(1 )rD 1 'j 

4 

2 

i-^.l 

(X*K 2 jIDj 

4 

2L X 


(1 + K 2) 

4 

2 

‘-^,1 

<I* K X,1> D 1 

4 

2Lj 

] - K X,X 

(1 + K x 1 )D 1 

4 

2Lj 

‘-*2,1 

(I +k 2,1> D X 

4 

2Lj 

1 _> 


(B19) 


For the thermal neutron flux 


<P 2 (X) = <p“(X) + <?£(X) 


= A 3 ex P^j + A 4 ex P ^ 


(B20) 


D l r l 


G 1 (X) = — L -± exp(-yX) - 

1 r\ - -rv 


D 2 n 


D„ 


A1 “ p @* A, “p(‘h)' 


x\ L l v l ex P(-^ X ) 


r 2 L 2 - X 


As in equation (Bll), 
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<Po (X) = -mxp(— 

\ L 2, 


/[/•' 


(X)exp^L^dX 


exp 


4a> 

\ L 2 y 


IdX 




s r l 5 


<P 2 (X) = A g exp 


(x\ + a 4 eX p^ — LLJMJJ 


V 1 


\ 2/ D 2 (r 2 L^ - l)exp(yX) 


V T 2r 2 

^ / rI j 1 Ij 2 


Aj exp 


© 


+ A 2 expf— 

\ L iy 


d 2 ( l ! - L f) 


(B21) 


(B22) 


where A^ is equal to A^ plus terms derived from the particular solution <P 2 (x). The 
constants Ag and A^ are determined by specifying the thermal neutron albedo at both 
ends of the reactor. 


•W - k 1,2 J -< X 1> 

(B23) 

J_< x 2> = K 2,2 J +< X 2> 

(B24) 

MX) d 2 dMx) 

(B25) 

J (X) - — - i A 

+ 4 2 dX 


<P 2 (X) D 2 dcpJX) 
J (X) - 2 + i 1 

(B26) 


4 + 2 dX 


The constants Ag and A^ are determined by simultaneous solution of these equations. 


A 


3 “ 


^10^12- Vll exp ( 

X 2 - X l\ 
L 2 / 

/ X 2 - X 1 
*0*11 exp ( 2 L 
- \ 2 > 

j - *7*12 



(B27) 
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where 


A 4 “ 


“ ^9 ^10 ex ^ 


^x 2 - x n 


V 9^11 exp (2 


x 2 - Xi > 


\ L 2 


- ^7^12 


exp 



(B28) 


^8 


^9 


^10 


!- K l,2 

< 1 + K 1,2> D 2 

4 

2L 2 

2 )«Pp< x 2 ) , 

(1 +K 2j 2 )D 2 d<p p (X 2 ) 

4 

2 ctX 

1 - K 2,2 + 

(1 + K 2 2 > D 2 

4 

2L 2 

>2 )<Pp(Xi) 

(1 +k 12 )d 2 d<p p (X 1 ) 

4 

2 dX 

_ 1-K 1,2 + 

(1 + Kj 2 )D 2 

4 

2 L 2 

_ 1-K 2 , 2 _ 

(1 + K 2|2 )D 2 

4 

2 L 2 


(B29) 


n L 2 


<Pp(X) = 


D i s r lJ 

■ T 

V 


y 2 Lj - l 

d 2 (>- 2l I - j ) 


exp(-yX) 


v t 2j 2 




A-^ expf— ) + Ag expf — 



d 2 ( l 1 - L ?) 


(B30) 
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d<P p {X) 


(B31) 


dX 



The fuel distribution required to obtain the desired power distribution G(X) is 


S f (X) = 


DjG(X) 


(B32) 


G(X) can be related to the exponential power distribution Y(x) derived in appendix A 
(eq. (A29)). From equation (B13) 


G(X) = yj exp(-yX) 


and from equation (B6) 


G(X) = 


vS f (X)<P 2 (X) 


Equation (B6) can be rewritten as 


G(X) = 


yS f <P 2 


Z f (X)cp 2 (X) 


s f(p 


r 2 


where Sj<p 2 is the axially averaged fission rate and 

\(X)cp 2 (X)\ 
S f^2 


(B33) 


is the local-to-average fission rate that is also equal to the local-to-average power 
density Y(x)/Y of equation (A29). 

The expression Y(x)/Y must be converted from the nondimens ional variable x to 
the dimensional form 
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X 


x = 


Therefore 


v S f <p ? 
G(X) = -L-±. 

Di 


P 


exp - 


1 X-Xi 


\ X - X 1 l\ 

^ x 2 - X 1 


V/ 3 X 2 " X l/ 


- exp 


r 1 


(B34) 


A comparison of equations (B13) and (B34) results in 


y 1 = vZ f <P 2 exp 


V = [«*2 - Xl>] 


-l 


(B35) 


(B36) 


The average fission rate can be expressed in terms of the radial component of the volu- 
metric heating rate H(r) as follows: 


Va- 



Equation (B36) becomes 



(B37) 
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APPENDIX C 


SYMBOLS 


A 1 

constant used in equations de- 
scribing fast and thermal neu- 

p 

tron flux, neutrons/(m )(sec) 

"P 

ratio of maximum to axially aver- 
aged power generation rate for 
reactor with "chopped cosine" 

a i 

coefficients in polynominal de- 


axial-power shape 


scribing /3 as function of tem- 
perature, K" 1 

P d 

total reactor power output with 
cosine axial-power distribution, 

c p 

specific heat of coolant, J/(kg)(K) 


W 

D 

neutron diffusion coefficient, m 

P e 

total reactor power output with 
exponential axial-power distri- 

De 

F 

equivalent diameter of coolant 
channel, m 

mass flow rate of coolant, kg/sec 

Pr 

Q 

bution, W 
Prandtl number 

volumetric heating rate expressed 

H(r) 

h 

radial power distribution in fuel 

3 

element, W/m 

forced convection film coefficient, 
W/(m 2 )(K) 

Q/A 

in dimensional space variables, 
W/m 3 

average heat flux from surface of 
fuel rod, W/m 2 

J + 

neutron partial current in positive 

p 

direction, neutrons/(m )(sec) 

q 

volumetric heating rate expressed 
in dimensionless space vari- 

J_ 

neutron partial current in nega- 


bles 


tive direction, 

2 

neutrons/(m )(sec) 

R 

R n 

dimensional radial field point, m 
outer radius of fuel element, m 

K iJ 

neutron albedo at position i for 


energy group j 

Re 

Reynolds number 

k 

thermal conductivity, W/(m)(K) 

r 

dimensionless radial field point 

L J 

neutron diffusion length for energy 
group j, m 

S 

wetted perimeter of coolant 
channel, m 

m(x) 

fraction of total fuel rod power 

T 

dimensionless temperature 


generated up to axial position x 

X 

dimensional axial field point 

Nu 

Nusselt number 

AX 

reactor length, m 

n(x) 

ratio of local to axially averaged 
power generation 

X 

dimensionless axial field point 
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I 

j 


Yj relative local power at inlet end 
of reactor 

Y(x) relative axial power distribution 

a j factor defined by eq. (A22) 

fi thermal-hydraulic factor defined 

by eq. (A13) 

y exponential coefficient used to de- 

scribe axial-power shape, 

e energy released per fission, J 

?? number of neutrons produced per 

thermal neutron absorbed in 
fuel 

9 dimensional temperature, K 

9-^ prescribed temperature limit for 

fuel, K 

v number of neutrons produced per 

thermal fission 

£ a macroscopic absorption cross 

section for thermal neutrons, 

-1 

m 

£-. macroscopic thermal- neutron 

1 _ 1 
fission cross section, m 


S D removal cross section for fast 

K -1 
neutrons, m 

<pj neutron flux for energy group j, 
neutrons/m^-sec 

ip arbitrary constant 

Subscripts: 

c coolant 

f fission 

M fuel 

R removal 

1 inlet end of reactor or fast- 

neutron energy group 

2 outlet end of reactor or thermal- 

neutron energy group 

3 extrapolated distance to zero neu- 

tron flux 

Superscripts: 

F fuel material 

N nonfuel materials 


3P 
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